Explosive percolation via control of the largest cluster 
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We show that only considering the largest cluster suffices to obtain a first-order percolation 
transition. As opposed to previous realizations of explosive percolation our models obtain Gaussian 
cluster distributions and compact clusters as one would expect at first-order transitions. We also 
discover that the cluster perimeters are fractal at the transition point, yielding a fractal dimension 
of 1.23 ± 0.03, close to that of watersheds. 
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Percolation, the paradigm for random connectivity, has 
since Hammersley ^ been one of the most often applied 
statistical models [11 [3]. Its phase transition being re- 
lated to magnetic models [4 is in all dimensions one of 
the most robust second-order transitions known. This ex- 
plains the enormous excitement generated by the recent 
work by Achlioptas, D'Souza, and Spencer [5 describ- 
ing a stochastic rule apparently yielding a discontinuous 
percolation transition on a fully connected graph. Subse- 
quent work applied the process on other networks [6HTT] . 
However, reported results of finite-size studies and size 
distributions are not consistent with a first-order transi- 
tion. Since then various rules have been devised p!2HT4] 
and even a Hamiltonian formalism was proposed [15], 
all attempting a discontinuous transition towards an in- 
finite cluster. In all proposed models one tries to keep 
the clusters of similar size and some authors additionally 
suppress the internal bonds of clusters [5l [15] . Could one 
obtain a clear and consistent first-order percolation tran- 
sition? It is the objective of the present Letter to answer 
this question. One criterion is the cluster size distribution 
at the percolation threshold. Radicchi and Fortunato [9] 
as well as Ziff [7 found a power-law distribution with 
an exponent close to two. Although, different from the 
exponent of classical percolation the sole fact of finding 
a power law is untypical for first-order transitions. Also 
unusual for a first-order transition is that the clusters are 
fractal, as we found happens for the Achlioptas rule, from 
the behavior of the order parameter with the system size 
[3 [9] . It is a purpose of the present Letter to present a 
model in which a Gaussian cluster size distribution and 
compact clusters can be achieved in a systematic way, 
characterized by a fractal perimeter yielding a fractal di- 
mension similar to the one of watersheds and random 
polymers in strongly disordered media. 

Usual bond percolation can be implemented on a 
square lattice by randomly occupying bonds between 
neighboring sites, reaching its threshold at a certain frac- 
tion when opposite borders are first connected through 
one large cluster [2l[T6l[l7]. This percolation threshold, 
is characterized by the continuous vanishing of the order 
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FIG. 1: (Color online) Snapshots of the system, obtained on 
a square lattice with 1024^ sites, at pc, for four different bond 
percolation models, namely, classical 2^, Achlioptas product 
rule [6], largest cluster model {a — 1), and Caussian model 
(a = 1). The largest cluster and Caussian models are intro- 
duced in this Letter. 



parameter, i.e., a second-order transition. On a fully con- 
nected graph, Achlioptas et al. [5 , used the best-of-two 
product rule studied in detail by Friedman and Lands- 
berg [12^. Ziff reported simulations on a regular square 
lattice [6] [7], while Radicchi and Fortunato [8l[9] and Cho 
et al. [W on scale-free networks. 

More recently, other approaches have been introduced 
to obtain explosive percolation. Instead of a best-of-two 
rule Manna et al [13 , Cho et al. [11 , and Moreira et 
al. [15 proposed a weighted rule where bonds are occu- 



pied according to a certain probability. However, despite 
being rejection- free schemes, they are hmited to smah- 
system sizes and/or reduced number of samples. Here, 
we suggest an acceptance method where new bonds are 
selected randomly and occupied according to a certain 
weight yielding, for the first time, a clear first-order tran- 
sition. The considered scheme allows to consider system 
sizes 64 times larger than before [13 , specifically, we con- 
sider systems of 4096^ sites and averages over 10^ sam- 
ples. 

In our simplest rule {^^ largest cluster model"), as for 
classical bond percolation, a link is randomly selected 
among the empty ones. If its occupation would not lead 
to the formation or growth of the largest cluster, it is al- 
ways occupied, otherwise, it is occupied with probability 



min < 1 , exp 



(1) 



where s is the size of the cluster that would be formed 
by occupying this bond and s the average cluster size af- 
ter occupying the bond. The parameter a controls the 
allowed size dispersion. Note that, for a < 0, since the 
size of the largest cluster is always greater (or equal) 
than the average cluster size, all new bonds are occupied 
reducing to classical bond percolation, characterized by 
a continuous transition at the percolation threshold [2 . 
For a > 0, the probability of Eq. ^ suppresses the for- 
mation of a cluster significantly larger than the average, 
inducing a homogenization of cluster sizes. The Gaus- 
sian function has been considered because this is what 
we expect for the cluster size distribution at a first-order 
transition. However, to observe a discontinuous transi- 
tion any other function could be chosen, as far as it con- 
strains the largest cluster differing significantly, in size, 
from the average cluster. 

For nonequilibrium problems, where a free energy can- 
not be defined, transitions can still be classified based on 
the behavior of the order parameter [18 . A first-order 
transition, is characterized by a jump in the order param- 
eter, otherwise, a transition is denoted as continuous. For 
percolation, we define as order parameter the fraction of 
sites in the largest cluster (Poo) [2 . Here we also consider 
two other quantities: the second moment of the cluster 
size distribution (x), defined as 



X- 



E^ 



(2) 



where the sum runs over all clusters z, and the standard 
deviation (xoo) of the largest cluster size (smax) over dif- 
ferent samples. 
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FIG. 2: (Color online) Size dependence, for the largest cluster 
model, of the susceptibility (x), fraction of sites in the largest 
cluster (Poo), and its standard deviation per site (xoo/A/") at 
the percolation threshold, on a square lattice of linear size 
(L) ranging from 32 to 4096. All bonds are occupied with 
the same probability except the ones that lead to the for- 
mation/growth of the largest cluster, to which an occupation 
probability q is assigned, Eq. (IT]), with a = 1. Results have 
been averaged over 10^ samples. 



To estimate the percolation threshold we consider the av- 
erage value of p (fraction of occupied bonds) at which a 
connected path linking opposite boundaries of the sys- 
tem is obtained. Gonsidering different system sizes, for 
a = 1, we obtain for the percolation threshold Pc = 
0.632 ± 0.002. To identify the order of the transition, 
in the largest cluster model. Fig. [2] presents a finite-size 
study for Poo, X? ^^d Xoo/^, averaged over 10^ samples 
of square lattices with linear sizes ranging from 32 to 
4096. As we can see in the top inset of Fig. [2J above a 
certain system size, the order parameter, at the percola- 
tion threshold, does not show any finite-size dependence, 
staying at a constant value in the thermodynamic limit 
(L -^ oo). The second moment of the cluster size dis- 
tribution (x) scales with L^{d = 2) which is a sign of a 
first-order transition [191 ED] • The standard deviation of 
the largest cluster (smax) per lattice site, which was also 
considered in Refs. [6] and [7], converges, for larger sys- 
tem sizes, to a constant value, corroborating the presence 
of a discontinuous transition. 

To explicitly control the cluster size distribution we 
also implemented the following model. A new bond is 
chosen from the list of empty ones and occupied with 
probability given by Eq. (IT]). For internal connections we 
consider s as twice the cluster size. Since equation (IT]) is a 
Gaussian with average size s and size dispersion s/y/2a^ 
we denote this model as Gaussian model. Note that here 
the occupation probability is assigned to all new bonds 
even when they are not related to the largest cluster. 
This not only guarantees the control over clusters greater 
than the average, as in the previous model, but also over 
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FIG. 3: (Color online) Size dependence, for the Gaussian 
model, with a = 1, of the susceptibility (x), fraction of sites 
in the largest cluster (Poo), and its standard deviation per 
site (Xoo/^) at the percolation threshold, on a square lattice 
of linear size (L) ranging from 32 to 4096. All bonds are 
occupied with a probability given by Eq. (fTl). Results have 
been averaged over lO'^ samples. 
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FIG. 4: (Color online) Number of sticks necessary to fol- 
low the perimeter of the infinite cluster as a function of the 
stick length, to obtain the fractal dimension of the perimeter 
with the yardstick method. For both the largest cluster and 
Gaussian models, with a = 1. For the Gaussian model data 
were vertically shifted by a factor of 0.1. Results have been 
averaged over 10^ samples of lattices with linear size 2048. 



the smaller ones. For a = 0, all bonds have the same 
probability and, therefore, the model reduces to classical 
bond percolation. For negative a, the growth of larger 
clusters is favored in two different ways: they differ more 
from the average value and have more empty bonds than 
the smaller ones. Yet, for all negative a, the model recov- 
ers the classical universality class of percolation [U [18] . 

As example, for positive a, we present, in Fig. [3J a 
size dependence study of the order parameter, second 
moment of the cluster size distribution, and standard de- 
viation per site of the largest cluster, for the Gaussian 
model, with a = 1, at the percolation threshold, on a 
regular square lattice with linear size (L) ranging from 
32 to 4096. Results were averaged over 10^ samples. We 
extrapolate, for the infinite system, a percolation thresh- 
old Pc = 0.56244 ± 0.00006. As for the largest cluster 
model, the density of the infinite cluster does not change 
significantly with the system size, the second moment of 
the cluster size distribution scales with L^{d = 2), and 
the standard deviation per site of the largest cluster con- 
verges to a non-zero constant. As before these results 
imply a first-order transition. 

Figure [T] shows snapshots for four different models of 
bond percolation: classical, product rule, largest cluster 
model, and Gaussian model. All figures have been ob- 
tained at their respective percolation thresholds (pc)- For 
classical percolation and for the product rule, clusters of 
very different sizes are obtained. In fact, the cluster size 
distribution is characterized by a power law [71 [9] . How- 
ever, for the largest cluster and the Gaussian model, a 
characteristic cluster size is observed. Both models lead 
to a localized cluster size distribution. Small size disper- 



sion and number of clusters are observed for the largest 
cluster model. According to Eq. ([l]), increasing the value 
of a decreases the size dispersion. 

As clearly seen in the snapshots of Fig. [l] clusters ob- 
tained with our models are compact but we find that the 
surface is fractal. For the Gaussian model, we calculate 
for the cluster perimeter a fractal dimension of 1.23±0.03, 
obtained with the yardstick method [21 (Fig. [4]). For the 
largest cluster model, it is also characterized by a fractal 
perimeter with a fractal dimension of 1.26 ±0.04 (Fig. [4]). 
Compact clusters with fractal surface were also reported 
for irreversible aggregation growth in the limit of high 
concentration by Kolb et al. [22 . For the present models, 
the percolation thresholds are larger than the ones from 
previous models due to the compactness of the clusters. 

In Fig. p] we see the cluster size distribution, P{s^a)^ 
for different system sizes, obtained with the Gaussian 
model. Measurements have been performed at the per- 
colation threshold on a square lattice with 1024^, 2048^, 
and 4096^ sites, and averaged over 10^ samples. Three 
characteristic peaks are observed. In fact, the third peak 
(around 0.7) is only due to the largest cluster and only 
appears due to the small number of clusters at the perco- 
lation threshold, being finite-size effect. This peak is not 
observed when we compute the same distribution without 
considering the largest cluster. In the thermodynamic 
limit, since an infinite number of clusters exists, the con- 
tribution of a single cluster to the distribution vanishes. 
The presence of two main peaks is characteristic for a 
first-order transition showing, for a finite system, at the 
percolation threshold, coexistence of the percolative and 
non-percolative states [23] . 
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FIG. 5: (Color online) Cluster size distribution for the Gaus- 
sian model for different system sizes (a = 1), at the percola- 
tion threshold, on a square lattice, averaged over 10^ samples. 
Black-dashed lines are two Gaussian distributions fitting the 
results from simulation. The black-solid line is the sum of 
both curves. 



In conclusion, the present work reveals that, to obtain 
explosive percolation on a regular lattice it is sufficient 
to control the formation and growth of the largest clus- 
ter, instead of applying a rule to the overall set of empty 
bonds. We propose the largest cluster model which sys- 
tematically suppresses the formation of a largest clus- 
ter. We introduce as well, the Gaussian model, where 
a weight is assigned to each selected bond, such that 
a Gaussian distribution of cluster sizes is obtained, re- 
vealing the coexistence of two states at the percolation 
threshold. Our models, yielding clear first-order transi- 
tions, show that explosive percolation can be obtained 
under less stringent conditions that previously thought 
shedding light on the minimum ingredients to trigger ex- 
plosive percolation. In fact, we believe that our restric- 
tions on the formation of a largest cluster differing sig- 
nificantly, in size, from the average, is the required nec- 
essary condition and hope that this statement can one 
day be formally proven. The value of the novel fractal 
dimension of percolation that we discovered in the clus- 
ter perimeters is intriguingly close to the one found for 
watersheds (1.211 ±0.001) [24^ and random polymers in 
strongly disordered media (1.22 ± 0.02) [25 . 
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